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Abstract 



We present a calculation of the single top quark cross section for proton- 
antiproton interactions with =1.8 TeV at the Fermilab Tevatron collider. 
We examine the effects of top mass, parton distribution functions, QCD scale, 
and collision energy, on each of the component production mechanisms, and 
study the kinematic distributions for standard model electroweak production. 
At the upgraded Tevatron with y/s = 2.0 TeV and high luminosity, it will be 
possible to test the nature of the Wtb coupling using single top production. 
We estimate the sensitivity to measure the single top cross section, and thus to 
directly measure Vtb and the top quark partial width. We show what happens 
to the Vtb measurement when an anomalous (V+A) component is added to 
the Wtb coupling, and how the top quark polarization affects the kinematic 
distributions. 
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INTRODUCTION 



Top quarks can be created via two independent production mechanisms in pp collisions. 
The primary mode, strong ti pair production from a gtt vertex, was used by the D0 and 
CDF collaborations to establish the existence of the top quark in March 1995 [[[],§. The 
second mode is electroweak production of a single top quark or antiquark from a Wtb vertex. 
This mechanism provides a sensitive probe for several standard model parameters. 

The experimental value of the top quark mass is 175 ± 6 GeV ||, which is in good 
agreement with the value of 177 ± 7 ijj GeV derived from electroweak measurements at 
LEP, SLC, SPS, Tevatron, and neutrino scattering experiments [f||. Since the mass is of the 
order of the electroweak symmetry breaking scale (vacuum expectation value = 246 GeV), 
it is a very promising place to look for deviations from the standard model [H . 

The coupling between a W boson and top and bottom quarks has not yet been studied 
directly. Top quark pair production is not the best process for probing this Wtb coupling 
or the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V t b since the top quarks 
are produced from a gtt vertex. Information about Wtb coupling from top quark decay 
is relatively inaccessible in ti pair production because the width of the decay into Wb is 
proportional to the branching fraction of t— >Wb, which is close to unity in the standard 
model (the top quark partial width Y{t^Wb) oc V t 2 b , where 0.9989 < |V t6 | < 0.9993 at the 
90% confidence level |7[]). Single top quarks are produced at hadron colliders mainly from 
a Wtb vertex, and thereby provide a direct probe of the nature of the Wtb coupling and of 
V tb . 

In hadron collisions, there are several partonic processes which produce single top quarks 
in the final state. The P^-gluon fusion mechanism q'g^-tqb has been studied in Refs 



The W* s-channel process q'q^tb has been examined in Refs. f21|-p7j for the Tevatron. Use 
of single top quark production as a tool for studying the Wtb coupling has been discussed in 
a number of papers for hadron colliders, including the Tevatron [p!T| , p^ , ^0| , ^ -|32l . Single top 



production has also been studied for e + e colliders p8|- |30| , |33| -|43|1 , and for &y interactions 

HhH. 

In this paper, we first present new results of consistent tree level cross section calculations 
for two and three vertex subprocesses of single top quark production in pp collisions at 
a/s = 1.8 TeV. We have chosen to use this energy since it is where the Tevatron collider 
operated between 1992 and 1996, and this work is in support of a search of the data for 
single top production. We then prepare the ground for single top physics at future high 
luminosity Tevatron runs with = 2.0 TeV, by making new estimates of the sensitivity to 
measure the top quark partial width, Wtb coupling, and V t b, including an anomalous (V+A) 
coupling. These studies are the first to be performed using complete tree level matrix element 
calculations for all possible processes, and do not use the effective W approximation method 

§■ 

In section 1, we provide a comprehensive overview of the three separate single top pro- 
cesses at the Tevatron, and their subprocesses. We describe our computation techniques, 
and present and discuss the results of the calculations. We have studied the cross section 
as a function of top quark mass, parton distribution parametrization, choice of scale Q 2 , 
and collider energy, and we have evaluated lower and upper bounds on the single top cross 
section. In section 2 of this paper, we investigate the kinematic distributions of single top 



2 



events, showing the separate contributions from the principal production modes. In section 
3, we look at the effects of a nonstandard coupling at the Wtb vertex, specifically the ad- 
dition of an anomalous right-handed (V+A) contribution to the coupling. We estimate the 
sensitivity to determine V t b by measuring the single top cross section at future Tevatron runs 
as a function of this possible right-handed coupling strength, and show how the polarization 
of the system may be used to help distinguish different scenarios. We also use our estimate 
of the cross section precision to predict the error on the top quark partial width at the 
upgraded Tevatron. Finally, in section 4, we summarize our results and draw conclusions. 



I. SINGLE TOP QUARK CROSS SECTION 

A. Single top processes 

Within the standard model, there are three separate processes at a proton-antiproton 
collider which result in a single top quark in the final state. The list below shows that 
these processes in turn consist of several subprocesses with two or three tree level vertices. 
The number of Feynman diagrams for each subprocess is shown. Some diagrams have been 
omitted from the total: (a) those with tt pair production from a gtt vertex and no electroweak 
vertex (not single top); (b) those containing a photon, Z boson or Higgs (their contribution 
to the total cross section is extremely small); and (c) diagrams with vertices containing 
off-diagonal CKM matrix elements. The notation used is that q is a light quark, and X 
represents any additional final state particles from the pp interaction. 

1. pp — > tb +X s-channel W* boson 

1.1 q'q-+tb_ 1 

1.2 q'g^tbq 2 

1.3 q'q^tbg 4 

2. pp — > tq + X t- or u-channel W boson 

2.1 q'b^tq 1 

2.2 q'g^tqb 2 

2.3 bg^tqq' 2 

2.4 q'b^tqg 4 

3. pp -> tW~ + X 

3.1 bg^tW _ 2 

3.2 qq-MWl 2 

3.3 gg^tWb 5 

3.4 bb^tWb 5 

3.5 qb-^tWq 3 

3.6 bg^tWg 8 

It should be noted that there is some variation in the literature over the use of the term 
"W-gluon fusion". In some papers it refers only to subprocess 2.2 q'g^tqb, in others to 
subprocesses 2.1 q'b^tq and 2.2 q'g^tqb combined. Only one reference flTj includes all four 



subprocesses 2.1-2.4 in the calculations, and the authors use the term 'W-gluon fusion" to 
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refer only to subprocess 2.2. In this paper, we will also use the term to mean only subprocess 
2.2. 

The two subprocesses 1.2 and 2.2 (both q'g^tqb), although superficially similar, each 
contain two different Feynman diagrams which are gauge invariant in separate pairs and do 
not need to be calculated together as a group of four diagrams. One might be tempted to 
consider the four diagrams as an independent set instead of as separate higher order cor- 
rections to the two main processes, since an experimental search will be able to distinguish 
between two-body and three-body final states. However, this is not acceptable mathemat- 
ically when calculating the cross sections, due to the definition of the b sea quarks in the 
parton distributions for the proton and antiproton. This definition requires the subprocesses 
to be grouped as shown in the list above. 

We will now discuss the subprocesses we have included in our calculations. In the fol- 
lowing list, the first initial state particle is a parton from the proton and the second one is a 
parton from the antiproton. Processes with an initial state c or s quark, or with off-diagonal 
CKM matrix elements are omitted from the list (and from plots) for simplicity, but have 
been included in our calculation of the overall cross section and other numerical results. 

1 . pp — > tb + X 

1.1 ud—>tb, du—>-tb 

1.2 ug—>-tbd, gu—>tbd, dg^-tbu, gd—>tbu 

2. pp — > tq + X 

2.1 ub—>td, bu—*td, db^tu, bd—>-tu 

2.2 ug—^tdb, gu^tdb, dg—^tub, gd-^tub 

3. pp -> tW~ + X 

3.1 bg^tW,_gb^tW 

3.2 uu^tWb, uu^tWb, dd^tWb, dd^tWb 

3.3 gg^tWb 

If the initial state parton is a u quark from the proton, then contributions from both valence 
and sea u quarks are included in the calculations. This also applies to d quarks in the 
proton and to antiquarks in the antiproton. Typical Feynman diagrams for these processes 
are shown in Fig. [I]. 

We have included in our calculations all the significant single top subprocesses with two 
or three vertices, except those with a gluon in the final state, which are significant, but 
which require a full next-to-leading order calculation to be included properly. Subprocesses 
with an extra quark in the final state, for instance bg^tqq' and qb—^tWq, although they 
have several Feynman diagrams, only contribute 1.5% to the total pp — > tq + X cross section 
and 1% to the pp — ► tW + X rate, and have therefore been ignored. This also applies to 
the subprocess bb^-tWb, despite its having multiple Feynman diagrams, including ones with 
electroweak ti production. 

B. Calculation details 

We have calculated the production cross section for each of the single top subprocesses 
mentioned in the previous section. We used the software package GompHEP B7[ to do the 
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tree level symbolic calculations and to generate optimized FORTRAN code for the squared 
matrix elements. Version 2.0 of CompHEP used the bases package to integrate over 
all phase space using parton distributions, and a CompHEP-BASES interface program to 



generate the correct event kinematics, with smoothing of singular variables [ 49[ . The Monte 



Carlo event generator SPRING j48| was used for each process in CompHEP 2.0. CompHEP 3.0 
has since replaced bases and SPRING with VEGAS | 50fl , and we used this version as well. 
Events generated were processed using pythia via a custom interface, in order to decay 
the W boson for use in kinematic studies of the final state particles. 

For these calculations, we have utilized the cteq3m and MRS (a') [^] parton distri- 
butions. These two sets of next-to-leading order structure functions both use the modified 
minimal subtraction (MS) renormalization scheme [54]]. The newly available parton distri- 
butions cteq4m |55| and mrs(r) are very similar to the distributions we have used, in 
the kinematic region for single top quark production at the Tevatron. 

We used the following standard model parameters in our calculations: Z boson mass 
rriz = 91.19 GeV, sin 2 6 W = 0.225, where 9 W is the weak mixing angle, (giving the W boson 
mass itl\y = mzCOs9 w = 80.28 GeV), b quark mass mj = 5.0 GeV, a = 1/128, and CKM 
matrix elements V u d = 0.975 and V t b = 0.999. All results have been obtained in two gauges, 
the unitarity gauge and the 't Hooft-Feynman gauge, as a check of calculations. Differences 
between calculations in the two gauges are less than 0.1%. 

We have chosen to use m 2 as the QCD evolution parameter or scale Q 2 value, since a 
large scale is a natural choice for top quark production. A high value such as m 2 is also a 
conservative choice that leads to lower cross sections. As shown in Refs. fOl and 111811, the 



leading order single top cross section depends rather strongly on the choice of QCD scale 
for values below ~(mn//2) 2 , and if a very small value for Q 2 such as m 2 were to have been 
chosen, then the calculated cross sections would have been about four times larger than 
those obtained using m 2 . 

A typical x value for single top quark processes is ~ mt/yfs ~ 180/1800 = 0.1, where x 
is the fraction of the proton or antiproton momentum carried by each initial state parton. 
At a scale Q 2 = (180 GeV) 2 , the value of a s is 0.102 from cteq3m and 0.104 from mrs(a'). 
Aqcd for five quark flavors is 158.0 MeV in both cteq3m and MRS (a'). 



C. Combining cross sections 

Care must be taken when combining some single top subprocesses in order to avoid 
double counting. One cannot simply add up the separate cross sections to get the total 
when there is a sea b quark in the initial state. The cteq3m and MRS (a') b sea distributions 
are not measured experimentally, but are obtained from the gluon distributions using the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations |p7 |. The b sea 



distributions in the structure functions therefore contain a mass singularity from the collinear 
divergence which occurs when the gluon splits to an onshell bb pair. The subprocesses we 
are considering which pertain to this situation are 2.1 q'b—>-tq, and 2.2 q'g^-tqb (H^-gluon 
fusion), where the initial state b quark in subprocess 2.1 is derived from the gluon sea in 
the antiproton. The correct way to avoid this singularity would be to calculate the rate for 
subprocess 2.2 q'g^-tqb with complete loop corrections, and then subprocess 2.1 with its b 
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sea contribution would be automatically included without the need for extracting it from the 
parton distribution sets. However, since we are making leading order calculations, we need 
another method. One technique for obtaining the cross section a(pp — > tq + X) is to 

calculate the rates for subprocesses 2.1 and 2.2 and add them together, and then subtract the 
rate from the splitting process g^bb convoluted with subprocess q'b—>tq. We have chosen 
to employ this method here. For instance, at m t = 180 GeV with Q 2 = m 2 and cteq3m 
parton distributions, the naive cross section for q'b—>tq is 0.75 pb, for q'g—+tqb it is 0.29 pb, 
and the splitting term is 0.54 pb, giving a total cross section from these two subprocesses of 
0.50 pb. The rate for pp — > tq + X is thus over 70% higher than the rate from the jy-gluon 
fusion subprocess q'g—^tqb alone. 

A second subtlety |Kj comes into play with this method for avoiding double counting, 
when working with the two parton distributions being considered: CTEQ3M and MRS (a'). 
The CTEQ collaboration has chosen to evolve the b sea distribution from m&, whereas the 
MRS group starts the evolution of its b sea distribution at 2m{,. The logarithmic terms 
in the gluon splitting must be evaluated as \n(Q 2 /m b ) with CTEQ distributions, but as 
ln(Q 2 /2rrib) for MRS distributions, to be consistent with the respective definitions of the 
b quark sea. 



D. Cross section versus top quark mass 



We show results for the cross sections of the three electroweak single top processes and 
the totals in Fig. || as a function of the top quark mass, with y/s = 1.8 TeV. Figure §(a) 
shows pp — > tb + X, (b) portrays the process pp — > tq + X, (c) is for the less important 
pp — > tW + X mode, and (d) shows the totals for each of these three processes for t and 
t combined. Figure 0(d) also shows ti pair production for comparison (upper line), and it 
can be seen that when only one top quark is produced in the final state, the cross sections 
decrease more slowly with increasing top quark mass than when two heavy tops have to 
be created at once. The strong ti cross section illustrated is from the resummed next-to- 
leading order calculation of Berger and Contopanagos |6lH , who used the cteq3m parton 
distributions. The tree level single top cross sections are the average of our calculations 
using cteq3m and MRS (a'). 

The main contribution to electroweak single top production comes from pp —>■ tq + X, 
the W boson t- and u-channel mode, including W^-gluon fusion. The rate from this process 
(61%, at m t = 180 GeV) is nearly twice as large as that from pp — > tb + X with a W* 
in the s-channel (32%). The contribution to the total cross section of the third process 
pp — > tW + X is small (7%). Of the dominant t- and u-channel process, 41% of the rate 
comes from q'b—^tq (after subtraction of the splitting term), and 59% from jy-gluon fusion 
q'g^tqb. Therefore, VT-gluon fusion forms 36% of the total single top rate from all processes. 

Because there are contributions from several single top processes, the total cross section 
forms a significant fraction of the ti pair production rate. The single top and antitop cross 
section from pp production at y/s = 1.8 TeV is 0.92 x 2 = 1.84 pb for a top quark of mass 
180 GeV and the cteq3m parton distributions, and 0.84 x 2 = 1.68 pb using MRS (a'). 
Therefore, although the rate of single top production is smaller than that from ti pair 
production (e.g. 4.7lig;35 pb at 180 GeV |61]| ) for all top quark masses considered here, it 
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is large enough to be extremely interesting for study at the Tevatron. 

Recent calculations show that higher order corrections to the leading order single top 
cross sections presented here are large. For instance, Ref. |17| shows that the K factor for 
pp — > tq + X is ~1.45 for m t = 180 GeV at y/s = 1.8 TeV with cteq2d parton distributions 
|62|j and Q 2 = m 2 . Reference [2^] contains a similar higher order calculation for pp — > tb+X, 



and finds that for m t = 175 GeV and y/s = 2.0 TeV, the K factor is also 1.45, using the 
CTEQ3M parton distributions and Q 2 = m^. 



E. Contributions to the single top cross section 

Table | presents values of various partonic subprocess cross sections for a top quark 
of mass 180 GeV, and for the two parton distributions discussed previously. Subprocesses 
with an initial state strange or charm sea quark contribute 1.9% to the total pp — >• tb + X 
cross section, and 6.1% to the total pp — > tq + X rate. Off-diagonal CKM matrix element 
subprocesses (not including initial state s and c sea quark subprocesses) contribute 0.3% 
to pp — > tb + X and 5.0% to pp —> tq + X. All these other modes contribute < 0.5% to 
pp — > tW + X production. Off-diagonal CKM subprocesses and initial state s and c sea 
quark subprocesses are included in our calculation of the total single top cross section, but 
are not included in our plots for simplicity, because of the calculation technique used. 

For single top modes like W* s-channel production with only light partons in the initial 
state, including both valence and sea quarks, the cross sections calculated with MRS (a') 
are 2.4% lower than those calculated with cteq3m. When there is a gluon in the initial 
state, for instance W-gluon fusion, then the mrs(a') cross sections are 5.7% lower than 
the cteq3m ones. For reactions with a sea b quark in the initial state, the cross sections 
calculated with MRS (a') are 17% smaller. 

Table |T| shows the resulting single top cross sections as a function of top quark mass. The 
central value numbers are the mean of the values calculated using cteq3m and MRS (a'). The 
upper and lower bounds come from combining half the difference between the calculations 
using the two parton distributions with the errors from the choice of QCD scale, as discussed 
in the next section. The correlation of the errors is correctly accounted for by adding the Q 2 
errors to each subprocess separately, before adding them in quadrature with the uncertainty 
due to the choice of structure function. For a top quark of mass 180 GeV, the total single 
top plus antitop cross section is 1.76i ;ig pb. 



F. Cross section versus scale 



We have examined the effect of the choice of QCD evolution parameter Q 2 on the various 
single top subprocesses. The results are shown in Fig. ^ for a top quark of mass 180 GeV and 
the cteq3m parton distributions at yfs = 1.8 TeV. Figure ||(a) shows the scale dependence 
for the W* s-channel process pp — > tb + X, which is the least dependent of the various 
single top processes on the choice of scale. Figure |3](b) is for the t- and u-channel processes 
q'g^tqb (W-gluon fusion) and q'b^-tq. The W-gluon fusion cross section falls rapidly as the 
calculation scale increases, whereas the subprocess q'b^tq goes up as Q 2 is raised. When 
these subprocesses are combined, the two effects partially cancel. The q'b^tq subprocess 
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is shown with the g^bb splitting term already subtracted. The minor single top process 
pp — > tW + X is shown in Fig. |3](c) with its various contributing subprocesses, which again 
have differing dependences on Q 2 that partially cancel in the sum. Finally, Fig. 0(d) shows 
each single top process summed for t and t, and the total single top production on the same 
axis for comparison. 

Leading order cross sections show more sensitivity to the choice of scale than higher order 
calculations. However, one still needs to choose a scale at which to perform the calculations. 
From an intuitive perspective, it does not make sense to consider top quark production 
as occuring at the almost zero mass b quark scale, although it would be mathematically 
consistent. Therefore, for our calculations, we have chosen the central value of the scale to be 
Q 2 = m 2 , and when estimating the uncertainty due to the choice of scale, we have restricted 
the region of interest to lie between (m t /2) 2 and (2m t ) 2 , as shown on the upper axes of Fig. 
0. The resulting errors on the main contributions to the cross section at m t = 180 GeV are: 
lxi% for pp — > tb + X; +9% for q'b^-tq and ^20% for q > g—>tqb, which combine to give i§ 5 % 
for pp — > tq + X; and tff% for pp — > tW + X. The q'b^tq and q'g^tqb errors largely cancel 
because the contributions to the errors from choice of Q 2 are 100% anticorrelated. The Q 2 
scale error dominates the total errors on the cross sections given in Table |H|. 

Combining the subprocesses q'b—>tq and q'g^tqb by subtracting the gluon splitting term 
to avoid double counting, as discussed earlier, is a procedure that is sensitive to the choice of 
evolution parameter Q 2 . Figure |] shows the subprocess q'b—*tq before and after subtraction, 
as a function of (a) the top quark mass, and (b) the scale Q 2 . It can be seen however, that 
provided the scale remains in the region around m 2 , then the sensitivity is less than that 
seen for the H^-gluon fusion subprocess q'g^tqb in Fig. 0(b). 



G. Cross section versus collider energy 

We have calculated the single top quark cross section as a function of production energy 
\/s. Figure 0(a) shows the cross section versus top quark mass for four collision energies: 
(i) the current Tevatron energy 1.8 TeV; (ii) the Tevatron energy for the next run in 1999, 
2.0 TeV; (iii) the energy of a possible Tevatron upgrade, 4.0 TeV; and (iv) the energy of 
the Large Hadron Collider (LHC) at CERN in 2005, 14 TeV. The three Tevatron cross 
sections are for pp collisions, whereas the LHC cross sections are calculated for pp collisions. 
Despite the ~150x increase in cross sections at the LHC, it will still be rather difficult to 
study single top quark production there, since the backgrounds will be much larger, and 
the signal will be harder to identify, because the jet produced at the same time as the top 
quark in VT-gluon fusion for instance, will be further forward in pseudorapidity rj, where 
i] = lntan(#/2) and 9 is the polar angle between the jet and the proton beamline. Peaks in 
the accompanying jet distribution at the LHC will occur at 77 = ±2.5 (c.f. 77 peaks at ±1.5 
when y/s = 1.8 TeV). 

At mt = 180 GeV, the cross section for single top quark production is 0.85 pb at 1.8 TeV, 
1.4 pb at 2.0 TeV, 9.4 pb at 4.0 TeV and for pp collisions at 14 TeV, 179 pb. For 180 GeV 
t antiquarks, the cross sections are the same as for t quarks at the Tevatron, but only 
133 pb at the LHC (26% lower), because there are no valence antiquarks in the initial state. 
These calculations were done using the CTEQ3M parton distributions with Q 2 = m|, and no 
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contributions from initial state s or c quarks or off-diagonal CKM matrix element terms are 
included. 

The relative contributions to the total single top cross section from each of the significant 
processes is not the same at all production energies. For a top quark of mass 180 GeV, Fig. 
|5](b) shows the single top plus antitop cross section versus production energy at the Tevatron 
for each component of the signal separately. It can be seen that the W* s-channel process 
pp —>■ tb + ib + X is much less sensitive to the change in available energy than the other 
processes, which increase rapidly in rate as the initial state energy goes up. At y/s = 1.8 TeV, 
the W* process forms 32% of the total single top signal, at 2.0 TeV it provides 29% of the 
cross section, and by 4.0 TeV it contributes only 13%. The pp — ► tb + ib + X process behaves 
in this manner because it is an s-channel process and its contribution to the total cross section 
comes from the s threshold phase space region, which is independent of energy. The reason 
why the total cross section increases with energy is that at higher energies, regions of smaller 
x in the proton structure functions are probed, and this is where the parton distributions 
are larger. 

The contribution from pp — > tW + X to the total single top cross section increases from 
7% at v/i = 1.8 TeV through 9% at 2.0 TeV, to 20% at 4.0 TeV. At the LHC, pp—>tW + X 
will contribute 30% of the single top rate and 40% of the antitop rate, and could therefore be 
an important production mode in the future. On the other hand, at the LHC, the s-channel 
process pp^tb + ib + Xwill fall to only 5% of the total single top rate, and will become 
experimentally inaccessible. 

H. A closer look at VF-gluon fusion 

We have analyzed the contributions to the production rate from the two Feynman di- 
agrams which form W^-gluon fusion, q'g—*tqb, shown in Fig. ||(a). There is no interference 
between the W^-gluon fusion diagrams and the two nonfusion diagrams of subprocess q'g^tbq 
(shown in Fig. |l](a) as subprocess 1.2), because the final state t quark and b antiquark have 
a different color structure. For the nonfusion diagrams, the t and b are from a W decay and 
so are in a color singlet state, whereas for the fusion diagrams the t and b come from a gluon 
and so are in a color octet state. 

The contribution to the total production rate of iy-gluon fusion from the Feynman 
diagram where the gluon produces a ti pair is very small, at about 5%. However, this 
diagram interferes destructively with the main VT-gluon fusion diagram where g-^bb. The 
destructive interference reduces the total rate for H^-gluon fusion by 34%. We present the 
cross section versus top quark mass for the two diagrams of H^-gluon fusion separately, and 
show the interference and net result, in Fig. §(b). 

I. More on pp —> tW + X 

We have considered two related 2^3 body processes in addition to the process bg^tW. 
These are qq^tWb and gg^tWb. We looked at these processes because in e + e~ and 77 
colliders, single top quark processes with iWb in the final state are important. However, 
we found that at the Tevatron these processes are not very significant. The interactions 
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e + e~, 77, qq, and gg^tWb all include diagrams with tt pair production and subsequent 
decay of the t into Wb, as well as many additional diagrams with just single top quark 
production. One needs to remove the contribution to the cross section from the invariant 
mass region m t = my/b around the top quark pole in order to study the Wtb vertex in single 
top production. The remaining contributions in e + e~ and 77 collisions are large enough 
(at 10 fb which is 17% of the total tWb cross section for e + e~ collisions at y/s = 2 TeV 
for example) to be sensitive to the coupling structure, but in qq and gg collisions almost 
the entire cross section comes from the tt diagrams, and the remaining single top quark 
contribution at 29 fb, is only ~0.8% of the total tWb rate of 3.5 pb. 

II. KINEMATIC DISTRIBUTIONS 

In order to understand in more detail the properties of single top quark production, 
we present in this section several experimentally interesting kinematic distributions. These 
are shown for top production only (not t) to make the presentation clear. Distributions 
for antitop are the same as those for top in transverse momentum, but are mirror images 
in pseudorapidity. If the sign of the W boson charge can be measured using its leptonic 
decay mode, then it will be possible to study the properties of top quarks and i antiquarks 
separately. All plots are for a top quark of mass 180 GeV and have been calculated using 
the cteq3m parton distributions at y/s — 1.8 TeV. 

The top quark decays to a W + boson and a b quark, and we consider here only subsequent 
leptonic decays of the W to a positron and neutrino, as this signal should be easier to find 
experimentally than channels with hadronic decay of the W boson. The branching fraction 
B for this decay mode is ~. The signature for a single top quark event is therefore a central, 
isolated, high px lepton), large missing transverse momentum from the neutrino, and at least 
two jets, where one of the jets comes from the hadronization of the b quark from the decay 
of the top quark. All single top events therefore have one potentially identifiable b jet, and 
of the experimentally accessible production modes at the Tevatron (pp — > tb + X , tq + X), 
~71% of them have a b jet as well. 

A. Transverse momentum 

Figure [j] shows the branching fraction times differential cross section B ■ da/dpx versus 
transverse momentum p T of the final state partons in single top production, and their decay 
products. In each plot, the short-dashed line is for W* production q'q^tb, the longer-dashed 
line for the two-body t-channel process q'b—>tq, and the narrow solid line for VT-gluon fusion 
q'g—>-tqb. The wide solid line is the sum of these three processes. Plot (a) shows the 
transverse momentum distributions of the top quark from each single top process. The 
mean of these distributions is 51 GeV. Despite its very high mass, the top quark is not 
produced at rest, but carries considerable transverse momentum in all three production 
modes. When the top decays, it produces a b quark, whose px distribution is shown in 
plot (b). The mean p T here is 62 GeV. Plot (c) is for the light quark produced with top 
in the t-channel processes ((px) = 43 GeV), and (d) is for the b antiquark often produced 
with top. Here the b from W* single top production has (px) = 59 GeV, whereas the b in 
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W^-gluon fusion is much softer, with (p?) = 25 GeV. The low will make this jet much 
more difficult to reconstruct. When the top quark decays, it produces a W boson, whose 
Pt is shown in Fig. 0(e) ((pt) = 65 GeV). The W decays to a positron, (shown in (f)) with 
mean p T of 45 GeV and a neutrino (in (g) , 48 GeV) . 



B. Pseudorapidity 

Figure |8] shows the branching fraction times differential cross section B ■ da/d-q ver- 
sus pseudorapidity r\ of the final state partons from single top production and their decay 
products. Plot (a) is for the top quark itself, where one can see that the pseudorapidity 
distributions are rather broad, and that the contributing production modes have very dif- 
ferent kinematics from each other. Both the VT-gluon fusion and W* modes produce top 
quarks more in the forward or +77 direction than backwards (with the distributions peaked 
at 77 ~ 1.7) whereas the two body t-channel process q'b^tq produces mainly backwards 
traveling top quarks, with the peak at 77 ~ —2.3. This distribution is also narrower than 
the other two. We see next that the decay products from the top are produced much more 
centrally. Plot (b) shows the b quark pseudorapidities. The distribution for the b from top 
decay in ly-gluon fusion is peaked at 77 ~ 0.1, and the b from top in W* production at ~ 0.2. 
The b from top in q'b^tq is still produced somewhat backwards, with a peak at ~ —0.8 
reflecting the direction of its parent. We would like to note that the 77 distribution of the 
b quark from the top decay in W^-gluon fusion is in agreement with that seen by C.-P. Yuan 
using the onetop generator but is rather different from the distribution for VT-gluon 
fusion shown in the TeV-2000 study of WH, H^bb |x| (with single top as a background), 
where the herwig generator |64| was used for this type of single top. Herwig seems to 



produce 6's in a symmetric peak in the region 1 < \rj\ < 5. This difference is not understood. 

One of the striking features of W-^-gluon fusion is the forward direction in which the light 
quark is produced [fTI]] . This can be seen in plot (c), where the light quark from q'g^tqb has 
a broad distribution, peaked at ~ 0.7. The effect is seen more emphatically in the two body 
t-channel mode where the peak occurs around 1.7, resulting in the summed distribution 
peaking at 77 ~ 1.5. The pseudorapidity distributions of the b antiquark produced together 
with top in 71% of single top events are shown in plot (d). Both distributions peak at 
77 ~ —0.4; the soft b from VT-gluon fusion has a rather broad spread in pseudorapidity, 
whereas the much harder b from W* production is produced in a narrower pseudorapidity 
peak. The r\ distributions of the W boson from the decay of the top quark, shown in plot 

(e) , are peaked at ~ 0.3 for W^-gluon fusion, at ~ 1.1 for W* production, and at ~ —1.2 for 
the q'b—>tq mode, echoing the directions of their respective parent top quarks. The positron 

(f ) and neutrino (g) distributions are more central versions of their parent W bosons. 



III. Wtb COUPLING AND V tb 

Since the top quark is rather heavy, we expect that new physics might be revealed at the 
scale of its mass. Many variants of nonstandard physics relating to this subject have been 
considered in the literature. Possible anomalous gluon-top quark couplings are discussed in 
Refs. p5|-[7I[]. Contact terms and new strong dynamics involving the top quark have been 
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studied in fl72"|-f?"7|l . The Ztt coupling will be inaccessible until a high energy e + e~ or ji + n~ 
collider is in operation. Studies of the Wtb coupling however, will be possible before then 
using single top production at the Tevatron. 

In this section we examine the effects on single top quark production and on its decay 
kinematics of a deviation in the Wtb coupling from the standard model structure, and 
we consider how this will affect a measurement of the CKM matrix element Vtb- In the 
standard model, the Wtb coupling is proportional to V t b and has (V-A) structure. As 
explained in the introduction, the cross section for single top quarks includes the Wtb 
coupling directly, in contrast to ti pair production. Therefore, single top production provides 
a unique opportunity to study the Wtb structure and to measure V t b- Experimental studies 
of this type are among the main goals of single top physics. Because high statistics will be 
required to make sensitive measurements, all the results given in the remaining subsections 
of this paper are for single top events produced in Runs 2 or 3 of the Tevatron; that is, from 
1999 onwards, with a collision energy of yfs = 2.0 TeV. 



A. Anomalous (V+.A) coupling 

As an example of a deviation from the standard model Wtb coupling, we introduce an 
additional contribution from a nonstandard (V+A) structure with an arbitrary parameter 
A r , where the subscript r refers to the right-handed current it represents. In the unitarity 
gauge, the Wtb coupling is given by: 

<hh [7,(l-75)+A r7 ,(l + 75 )] 



2 v / 2sin# 

where e is the positron electric charge, sin6 w = 0.474, and 7 M and 75 are Dirac matrices. 

The dependence of the total single top quark cross section on the parameter A r is shown 
in Fig. |, for yi = 2.0 TeV, m t = 180 GeV, and V tb = 0.999. Here, a(pp -> t + i+ X) = 
2 x cr(pp — > tb + tq + tqb). The standard model value of A r is zero. The production rate 
varies almost quadratically with A r , and is nearly symmetric about the point A r = 0. The 
cross section rises from 2.44 pb when A r = to 4.68 pb when A r = —1 and to 4.73 pb when 
A r = +1. 



B. Sensitivity in the (V t b,A r ) plane 

We have calculated the region in the (V t b,A r ) plane for which there will be experimental 
sensitivity using future single top measurements. If one finds a number of single top events 
consistent with the standard model prediction, then it may be that the Wtb coupling is 
purely left-handed, and that V t b is close to unity. Alternatively, the cross section could be 
boosted by an anomalous contribution to the Wtb coupling, as shown for example in Fig. |9|, 
with Vtb correspondingly lower. 

The error on the measurement of Vtb is dependent on the error on the single top cross 
section, including both experimental and theoretical contributions. First we estimate the 
experimental error for a top quark of mass 180 GeV at y/s = 2.0 TeV as follows: we take the 
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integrated luminosity for Tevatron Run 2 as 2 fb 1 , with an error of 5%; the signal acceptance 



including at least one b tag as 0.20, from the TeV-2000 study of single top production |78 



with an error of 7%; and a signal to background ratio of 1:2, with a systematic error on the 
background of 7%. The available branching fraction includes both the electron and muon 
decay channels, giving B — |. In Run 2, all accessible modes of single top production will 
have to be used together in order not to make a statistics- limited measurement. We use here 
the value 2.44 pb for the single top cross section (= 0.72 pb (s-channel) + 1.72 pb (t- and 
u-channels)) from cteq3m. These assumptions lead to a prediction that approximately 650 
events will be found in a search, with one third coming from single top production and two 
thirds from various backgrounds (e.g. W + bb, W+ light jets with a mistag, tt). Therefore, 
the experimental error on the total single top cross section will be 10% (statistical) © 16% 
(systematic) = 19%, where the © symbol means "add in quadrature". 

The error on the theoretical calculation of the cross section includes contributions from 
the choice of parton distribution function and from the scale, as discussed earlier in this 
paper, where they were found to be ~12%. However, there is another contribution, not 
well quantified, from the lack of knowledge of the gluon distribution in the proton and 
antiproton for t- and u-channel single top processes. This error has been variously reported 
to the authors as 30% J7S| and 10% |81J, and so we use these values here to estimate the 
error on the theoretical total single top cross section at 32% or 16%. 

The error on a measurement of V t b will be half the error on the single top cross section, 
since the cross sections for all single top processes are proportional to \V t b\ 2 - This results 
in an error on V tb of (19% © 32%)/2 = 19% or (19% © 16%)/2) = 12% from the Tevatron 
Run 2, depending on one's view of the knowledge of the gluon momentum distributions in 
the proton sea. 

There may be a Run 3 at the Tevatron from 2002 onwards, producing 30 fb" 1 of data. 
This high luminosity mode of collider running is known as "TeV33" after the planned in- 
stantaneous luminosity of 10 33 cm _2 s _1 . With such high statistics available, it has been 
shown by Stelzer and Willenbrock [22] that using just s-channel W* production with double 
b tagging instead of all single top modes with only one tag will eliminate most of the un- 
certainty on the theoretical cross section, because there will no longer be any contributions 
from processes with initial state gluons. They also showed that the measurement should 
be possible using 3 fb -1 of data. We update their calculation here for Run 3, including 
estimates of the systematic errors. The cross section for pp — > tb + ib + X is 0.716 pb, with 
m t = 180 GeV, y^i = 2.0 TeV and Q 2 = m 2 . To estimate the error on Vtb using Run 3 
data, we make the following assumptions: the error on the luminosity remains at 5%; the 
signal acceptance for W* single top is 0.08 when requiring a double b tag, as shown in Ref. 



f22|1 , with a 1.8% error; and the signal to background ratio is 1:2 (again from [22]|), with 
a systematic error on the background of 1.8%. Therefore, an experiment at the Tevatron 
in Run 3 will see approximately 1,146 events when searching for W* single top production, 
with one third signal events and two thirds coming from various backgrounds (e.g. W + bb, 
W+ light jets with two mistags, WZ with Z^-bb, H^-gluon fusion, tt). This observation will 
lead to a measurement of the W* single top production cross section with an error of 7% 
(statistical) © 6% (systematic) = 10%. Smith and Willenbrock JM] show that the error on 



the theoretical cross section for W* single top production is only 3%, leading to an error on 
V lb of 
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In Fig. [K] we show the results of these calculations, extended into the (Vtb,A r ) plane. In 
plot (a) for Tevatron Run 2 (2 fb -1 ), the outer short-dashed contours show the result when 
the error on the theory cross section includes a 30% contribution from lack of knowledge 
of the gluon distribution. The inner long-dashed contours result from when this error con- 
tributes only 10% to the overall measurement. Plot (b) presents our estimates for Tevatron 
Run 3, "TeV33" (30 fb -1 ), with the dashed contours showing the precision obtainable using 
a theory error of only 3% and an experimental search to isolate the W* s-channel mode 
of single top. We discuss in the next section how one might distinguish standard model 
production from the scenario discussed above where the effects of an elevated cross 

section caused by the anomalous coupling cancel with a reduced value of V t b from a possible 
mixing of the top quark with a new fourth generation quark to give an observed number of 
events consistent with the standard model. 

C. Polarization of the top quark 

Top quark polarization depends strongly on the structure of the Wtb coupling, and one 
might expect an asymmetry in angular distributions of the final state partons for different 
values of A r . For example, standard model single top is produced almost 100% left-handedly 
polarized because of the left-handed current structure of the Wtb coupling, whereas if A r = 1, 
the top quark is not polarized at all. To calculate polarization effects using Monte Carlo 
generators, it is necessary either to keep the polarization of all particles in the final states of 
the 2— >2 and 2— >3 processes being studied (i.e. q'q—>tb, q'b—>tq, q'g—>tqb), with subsequent 
decays of the polarized top quark and W boson (t— >Wb, W^eu), or else one needs to 
calculate the higher order 2— >4 and 2— >5 processes (i.e. q'q^eubb, q'b^evbq, q'g^eubqb), 
with the top quark and W boson treated as resonances in the intermediate states. The 
second method automatically includes the polarizations of the intermediate state t and W. 
To study the differences between kinematic distributions when the polarizations of the top 
quark and W boson have been taken into account with those where they are assumed to be 
unpolarized (as in most Monte Carlo generators, e.g. pythia), we have calculated the 2— >4 
and 2— >5 processes for the three significant single top production modes using CompHEP 
alone, and compared the results with calculations where we used CompHEP for the 2— >2 and 
2^3 single top processes, and pythia for the subsequent t and W decays. 

Our direct calculations show that the px and rj distributions are not sensitive to the 
polarization of the top quark. 

Two representative examples of distributions expected to reflect the top quark polariza- 
tion effects are the invariant mass of the positron and the b quark, m e b, and the cosine of 
the polar angle, cos 9*. The invariant mass m e b is given by: 

m eb = yj (E e + E b ) 2 - (p Te + p Tb ) 2 - (p ze + p zb ) 2 

where p z is the momentum of the positron or b quark along the beam direction. 

The polar angle 9* is defined as the angle between the positron direction and the x axis 
within the rest frame of the W boson, where the x axis is defined to be in the direction of 
motion of the W boson in the rest frame of the top quark |16| . The cosine of this angle is 
given approximately by: 
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Figure |TT| shows the distributions of (a) m e & and (b) cos6**, for the case when polarizations 
of the top quark and W boson have been properly taken into account (solid histogram), and 
for when summation over the polarization of the top quark decay products has been done 
using the subsequent decays of an unpolarized top quark and W boson (dashed histogram). 
One can see that there are indeed differences in these distributions for the polarized and 
unpolarized cases. In particular, an asymmetry (or lack of it) in cos#* when the positron is 
emitted aligned or antialigned with the direction of motion of the W boson in the top rest 
frame should be observable with high statistics. All three modes of single top production 
exhibit this same behavior. The two variables m e b and cos 6* can also be used in combination 
with the total single top production rate, which is sensitive to the Wtb coupling structure 
as shown previously, to further our understanding of the Wtb coupling. 



D. Top quark partial width 

From our previous estimates of the sensitivity for measuring the single top cross section 
at future Tevatron runs, we can obtain the expected precision on the top quark partial width 
r(t— >Wb). In the standard model, the top quark partial width is very nearly the same as 
its full width, since V t b is so close to 1. The top partial width is directly proportional to the 
single top cross section, and so the error on the width is just the experimental measurement 
error on the cross section added in quadrature with the theoretical calculation error. Thus 
in Run 2 the top quark width should be measured with an error of 19% © 32% = 37% or 
19% ©16% = 25%, depending on whether the uncertainty on the gluon distribution function 
is 30% or 10%. In Run 3, this precision can be improved to 10% © 3% = 10%, which is 



comparable to what can be achieved at a linear e + e collider [31] using a tt threshold scan 



IV. CONCLUSIONS 

In this paper we have reported the results of new studies of single top quark physics 
at the Fermilab Tevatron pp collider. We have made consistent calculations of the tree 
level cross sections for each mode of single top production as a function of top quark mass, 
parton distribution function, QCD scale, and collision energy. We discussed details of the 
calculations for several of the subprocesses involved, and gave breakdowns of the various 
contributions to the overall cross sections. For a top quark of mass 180 GeV, at ^/s = 
1.8 TeV, with Q 2 = mf, and taking the mean result from cteq3m and MRS (a'), we find 
that the leading order total single top plus antitop cross section is 1.76i$jf pb. 

We have shown for each subprocess separately the transverse momentum and pseudo- 
rapidity distributions of the top quark, the other quarks produced with it, and its decay 
products. These kinematic distributions need to be understood in order to be able to sepa- 
rate signal from background in an experimental search. 

We then considered the possibility for measuring the CKM matrix element Vtb and the 
Wtb coupling directly using single top events from the next Tevatron run. We estimated 
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the sensitivity such measurements might have, and how an anomalous (V+A) term in the 
Wtb coupling would affect the measurement of V t b- If there is no anomalous component 
to the Wtb coupling, then V th can be measured to a precision of 19% or 12% in Run 2 
(1999-2001), with the two values coming from different estimates of the uncertainty in the 
gluon distribution function. In Run 3 (2002-2006), the precision on V t b will be improved 
to 5%. The top quark polarization affects the angular distributions of its decay products, 
and we investigated how this could be used together with a measurement of the single top 
cross section to distinguish between various processes affecting the top quark beyond the 
standard model. Finally our estimates of the single top cross section error show that the 
top quark partial width will be measured to within 37%-25% in Run 2, and to a precision 
of 10% in Run 3. 

We find the prospects for single top physics at the Tevatron exciting and that a rich 
program of studies will be possible in the future. 
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TABLE I. Production cross sections for sin- 
gle top quark processes for nit = 180 GeV, with 
Q 2 = m 2 and sfs = 1.8 TeV, using cteq3m and 
MRS (a') parton distributions. "Other modes" 
refers to subprocesses with an s or c quark in 
the initial state, and to subprocesses involving 
an off-diagonal CKM matrix element term. 
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TABLE II. Single top plus antitop cross sections at the Tevatron with yfs = 1.8 TeV and 
Q 2 = mf, as a function of top quark mass. The central values are the mean of the calculations 
using cteq3m and MRS (a'). The upper and lower bounds include the effects from choice of scale 
Q 2 and half the difference between the parton distribution functions. 
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FIGURES 




(c) 3.1 bg -> tW 



FIG. 1. Representative Feynman diagrams for the three single top quark production processes 
at the Fermilab Tevatron: (a) the W* s-channel process pp — > tb + X; (b) the W t- and u-channel 
process pp — > + X, including subprocess 2.2, IF-gluon fusion; and (c) pp — > tIF - + X. 
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FIG. 2. Single top quark cross sections at the Tevatron with ^/s = 1.8 TeV, versus top quark 
mass: (a) s-channel W* production pp — > tb + tbq; (b) t- and u-channel production pp — > tq + tqb; 
(c) pp — > tW + tWb; and (d) the total single top and antitop cross section pp — > t + t + X. The 
resummed next-to- leading order ti cross section of Ref. |6l| is shown as the uppermost line in (d), 
for comparison with single top production (at leading order). 
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Scale Q 2 [GeV 2 ] 

FIG. 3. Single top quark cross sections (mt = 180 GeV, \/i = 1.8 TeV) versus QCD evolu- 
tion scale Q 2 for: (a) s-channel W* production pp — ► tb + i&g; (b) t- and u-channel production 
pp ^ tq + tqb; (c) — > tW + iW6; and (d) the summed single top and antitop cross section 
pp-^t + t + X. 
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Top Quark Mass m t [GeV] Scale Q 2 [GeV 2 ] 



FIG. 4. Single top produced together with a light quark, q'b^tq, from an initial state sea b 
quark, showing the cross section before and after subtraction of the gluon splitting term, as a 
function of: (a) top quark mass (with Q 2 = mf); and (b) scale Q 2 (with m t = 180 GeV). 




FIG. 5. Single top quark cross section plotted (a) versus top quark mass, at four production 
energies: the Fermilab Tevatron at ^/s = 1.8 TeV; the upgraded Tevatron at 2.0 TeV; the proposed 
TeV* collider at 4.0 TeV; and the CERN pp Large Hadron Collider at 14 TeV. Plot (b) shows the 
cross section versus collider energy (with mt = 180 GeV), for each of the single top production 
mechanisms. The values in (b) up to 12 TeV are for pp production, whereas the results at 14 TeV 
are for pp collisions. 
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FIG. 6. (a) Feynman diagrams for P^-gluon fusion (q 1 g^tqb). (b) W-gluon fusion cross section 
versus top quark mass, showing the contributions from each of the Feynman diagrams, and the 
large destructive interference between the two processes. 
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FIG. 7. Single top transverse momentum distributions (with mt = 180 GeV) for: (a) the top 
quark; (b) the b quark from the decay of the top; (c) the light q quark produced with top in t- and 
u-channel processes; (d) the b antiquark produced with top in the s-channel W* process, and in 
T^-gluon fusion; (e) the W boson from the top decay; (f) the positron from the W decay; and (g) 
the neutrino also from the decay of the W boson. 



FIG. 8. Single top pseudorapidity distributions (with m t = 180 GeV) for: (a) the top quark; 
(b) the b quark from the decay of the top; (c) the light q quark produced with top in t- and 
u-channel processes; (d) the b antiquark produced with top in the s-channel W* process, and in 
VF-gluon fusion; (e) the W boson from the top decay; (f) the positron from the W decay; and (g) 
the neutrino also from the decay of the W boson. 
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FIG. 9. Total single top and antitop production cross section at the upgraded Tevatron, with 
yfs = 2.0 TeV and nit = 180 GeV, versus the right-handed (V+A) coupling parameter A r . 
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FIG. 10. Estimated la measurements in the (Vtb,A r ) plane for an experiment running at the 
upgraded Tevatron collider at \/s = 2.0 TeV, assuming that the number of events seen is consistent 
with the standard model prediction. Plot (a) shows the results for 2 fb _1 of data, using all accessible 
modes of single top production (pp — > tb + tq + tqb+ c.c). The outer short-dashed lines enclose the 
region resulting from a 32% error on the theoretical cross section and the inner long-dashed lines 
from a 16% uncertainty. Plot (b) is for a future run with 30 fb™ 1 of data, using only W* single top 
production (pp — * tb + ib + X) where the error on the theoretical cross section is 3%. 
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FIG. 11. Single top distributions of (a) invariant mass m e b, and (b) cosine of the polar angle 
6* (denned in the text). The solid histograms are for the standard model case where the top quark 
and W boson are ~100% left-handedly polarized (fully calculated using CompHEP for the 2^4 
and 2— >5 processes with intermediate state t and W resonances), and the dashed histograms are 
for when there is no polarization, corresponding either to a (V+.A) term with A r = 1 in the Wtb 
coupling, or to the case where the polarization has been excluded from the calculation (e.g. by 
using pythia to decay the W boson). 
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